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Abstract. Observations made using the Rossi X-ray Timing Explorer 
r^' have shown that accreting weak-field neutron stars in low-mass X-ray 

, binary systems (LMXBs) produce three distinct types of millisecond X- 

ray oscillations that can be used to determine the spin rates and estimate 
the magnetic fields of these stars. These oscillations show that more than 
two dozen neutron stars in LMXBs have spin rates and magnetic fields in 
the range that will make them radio-emitting millisecond pulsars when 
^> ' accretion ceases, supporting the hypothesis that neutron stars in LMXBs 

. are the progenitors of the rotation-powered millisecond pulsars. The spins 

of the 16 known accreting millisecond pulsars in LMXBs are consistent 
QQ . with spin-up to accretion spin equilibrium for magnetic fields ranging 

\ from 3 X 10^ G to 3 x 10^ G and time-averaged accretion rates ranging 

'nI" . from 3 x 10^'^ M^; to Mg, provided these stars have been accreting long 

enough. The P~P distribution of rotation-powered millisecond pulsars 

■ indicates that their initial spins are set by spin-up to spin equilibrium at 
accretion rates ranging from IQ-^Me to Me or that they never reached 

Q ■ accretion spin equilibrium. 

> ■ 

•'^ ■ 1. Introduction 

>< 

■ Soon after the discovery of millisecond pulsars (Backer et al. 1982), it was sug- 
gested (Alpar et al. 1982; Radhakrishnan and Srinivasan 1982) that they have 
been recycled by being spun up by accretion of angular momentum in low-mass 
X-ray binary systems (LMXBs). Although this picture has been widely accepted 
for more than two decades (see Bhattacharya 1995), until the last few years there 
was only indirect evidence that neutron stars in LMXBs have the spin rates and 
magnetic fields required for them to be the progenitors of the rotation-powered 
millisecond pulsars (MSPs). The Rossi X-ray Timing Explorer (RXTE) has 
transformed the situation by revealing that accreting neutron stars in LMXBs 
produce three distinct types of millisecond X-ray oscillations that can be used 
to determine their spins and estimate their magnetic fields. All three types of 
oscillation are thought to be generated directly or indirectly by the star's mag- 
netic field and spin. Together, they provide strong evidence that more than two 
dozen neutron stars in LMXBs have the millisecond spin periods and ~ 10^- 
10^0 G magnetic fields necessary for them to become rotation-powered MSPs 
when accretion ceases. 
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Periodic accretion-powered X-ray oscillations with millisecond periods have 

been detected in five neutron stars in LMXBs (see Tabic 1), establishing beyond 
any doubt that these stars have dynamically important magnetic fields. Their 
coherent oscillations indicate that they have field strengths ^ 10^ G (see Miller, 
Lamb, &; Psaltis 1998, hereafter MLP98), while their nearly sinusoidal waveforms 
and unusually low oscillation amplitudes indicate that they have field strengths 
^ 10^^ G (Psaltis & Chakrabarty 1999). The spin frequencies of these accretion- 
powered MSPs range from 185 Hz to 435 Hz. 

Nearly periodic X-ray oscillations have been detected during thermonuclear 
bursts of 13 neutron stars in LMXBs, including 2 of the 5 known accretion- 
powered MSPs (Table 1). The existence of thermonuclear bursts indicates field 
strengths ^ 10^'^ G (Lewin, van Paradijs, & Taam 1995) while the spectra of 
the persistent X-ray emission (Psaltis &; Lamb 1998) and the properties of the 
burst oscillations (see § 2) indicate field strengths ^ l(f G. The spin frequencies 
of these nuclear-powered MSPs range from 270 Hz to 619 Hz. 

Pairs of accretion-powered kilohertz quasi-periodic X-ray oscillations have 
been detected in more than two dozen accreting neutron stars (see Lamb 2003) , 
including 8 nuclear-powered MSPs and 2 of the 5 known accretion-powered MSPs 
(Table 1). The frequencies of these kilohertz QPOs range from ~ 100 Hz to 
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^ 1300 Hz, showing that accreting gas orbits close to the surfaces of these stars 
and that they have magnetic fields ~ 10^-10^ G (MLP98). The separation of 
the two QPOs remains constant to within a few tens of Hz as their frequencies 
vary by as much as a factor of 5 and is commensurate with the spin frequency of 
the star (see § 2), showing that the star's spin plays a central role in generating 
the QPO pair and indicating that the magnetic fields of these stars are ^ 10^ G. 
The spin frequencies of the stars that produce kilohertz QPOs can be inferred 
directly from the frequency of their periodic X-ray oscillations, if they have 
been detected, or indirectly from the separation of their kilohertz QPOs, if they 
have not. The spins of the neutron stars in which periodic X-ray oscillations and 
kilohertz QPOs have both been detected range from 191 Hz to 619 Hz (Table 1). 

These discoveries have established that many neutron stars in LMXBs have 
magnetic fields and spin rates similar to those of the rotation-powered MSPs. 
The similarity of these stars to rotation-powered MSPs strongly supports the 
hypothesis that they are the progenitors of rotation-powered MSPs. After being 
spun down by rotation-powered emission, the neutron stars in these systems are 
spun up to millisecond periods by accretion of matter from their binary com- 
panions, eventually becoming nuclear- and accretion-powered MSPs and then, 
when accretion ends, rotation-powered MSPs. 

In § 2 we describe in more detail the new evidence that many neutron stars 
in LMXBs have millisecond spin periods and dynamically important magnetic 
fields. In § 3 we discuss the spin evolution of weak-field neutron stars in LMXBs, 
including the long-standing idea that their spin rates are limited by coupling of 
their magnetic fields to the surrounding accretion disk and recent proposals that 
their spin rates may be affected by gravitational radiation by the spinning star. 
In §4 we summarize the implications of the new evidence for the evolution of 
neutron stars in LMXBs and formation of rotation-powered MSPs via recycling. 

2. Measuring the Spin Rates of Accreting Weak-Field Neutron Stars 

Neutron stars in LMXBs are accreting gas from a Keplerian disk fed by a low- 
mass companion star. The star's magnetic field and accretion rate are thought 
to be the most important factors that determine the accretion flow pattern 
near it and the spectral and temporal characteristics of its X-ray emission (see 
MLP98). The accretion rates of these stars vary with time and can range from 
the Eddington critical rate Me to less than 10~^Me- Their magnetic fields are 
thought to range from 10^^ G down to 10^ G or possibly less. Magnetic fields at 
the upper end of this range are strong enough to terminate the Keplerian disk 
well above the stellar surface, even for accretion rates ~ Me, whereas magnetic 
fields at the lower end of this range affect the flow only close to the star, even 
for accretion rates as low as ~ W~'^Me- 

For intermediate field strengths and accretion rates, some of the accreting 
gas is expected to couple to the star's magnetic field well above the stellar surface 
and be funncled toward the magnetic poles, heating the stellar surface unevenly. 
The remainder of the accreting gas is expected to remain in a geometrically thin 
Keplerian flow that penetrates close to the stellar surface, as shown in Figure 1. 
This flow is thought to be responsible for generating the kilohertz QPOs (see 
Lamb & Miller 2001; Lamb 2003; Lamb & Miller 2004). Recent observations (see 
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Figure 1. Side view of a weak- field neutron star accreting from a disk, 
showing the complex flow pattern expected. Some accreting gas couples 
strongly to the magnetic field and is funneled toward the magnetic 
poles, but a substantial fraction couples only weakly and drifts inward 
in nearly circular orbits as it transfers its angular momentum to the 
star via the stellar magnetic field. Prom MLP98. 



below) show that X-ray bursts also heat the stellar surface unevenly, and that 
the regions heated in this way are coupled to the rotation of the star, probably 
via the star's magnetic field. Whether due to accretion or to nuclear burning, 
uneven heating of the stellar surface produces a broad pattern of X-ray emission. 
Rotation of this pattern makes both the accretion-powered and nuclear-powered 
X-ray emission of the star appear to oscillate at its spin frequency. 

A breakthrough was achieved during the past year with the detection of 
burst oscillations in the accretion-powered MSPs SAX J1808.4— 3658 (Chak- 
rabarty et al. 2003) and XTE J1814-338 (Strohmayer et al. 2003) and precise 
measurements of the frequencies and phases of these oscillations. The new results 
show that, except during the first seconds of some bursts, fhe burst oscillations 
of these stars have the same frequency and phase as their coherent acc:retion- 
powered oscillations (see Fig. 2), establishing beyond any doubt that (1) these 
stars have magnetic fields strong enough to channel the accretion flow and en- 
force corotation of the gas heated by nuclear burning and (2) the nuclear- and 
accretion-powered oscillations are both produced by spin modulation of the X- 
ray flux from the stellar surface. The burst oscillations of some other stars are 
very stable (Strohmayer &; Markwardt 2002), but many show frequency drifts 
and phase jitter (Strohmayer et al. 1996; Strohmayer et al. 1998; Muno, Fox, 
& Morgan 2000; Muno et al. 2002). The new results confirm that burst and 
persistent oscillations both reveal directly the spin frequency of the star. 

The 16 MSP spins measured to date are consistent with a uniform distribu- 
tion that ends abruptly at 760 Hz (Chakrabarty et al. 2003), but they are also 
consistent with a distribution that decreases more gradually. The proportion 
of accretion- and nuclear-powered MSPs with frequencies higher than 500 Hz is 
greater than the proportion of known rotation-powered MSPs with such high fre- 
quencies, probably because there is no bias against detecting accretion-powered 
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Figure 2. An X-ray burst from SAX J1808.4-3658 observed on 18 
October 2002. Dark curve and scale at right: X-ray count rate as a 
function of time during the burst. Contours and scale at left: Dynamic 
power spectrum of the X-ray brightness, showing the rapid increase in 
the frequency of the burst oscihation from 397 Hz to 403 Hz during the 
rise of the burst, the disappearance of the oscihation at the peak of the 
burst, and its reappearance about 10 s after the start of the burst. The 
horizontal dashed hne shows the frequency of the neutron star's spin 
inferred from its accretion-powered periodic X-ray osciUations. From 
Chakrabarty et al. (2003). 



MSPs with high frequencies, whereas detection of rotation-powered MSPs with 
high frequencies is stiU difficult (see Chakrabarty et al. 2003). 

Seven of the 16 known MSPs in LMXBs have frequencies > 435 Hz, whereas 
none of the accretion-powered MSPs have frequencies this high (see Table 1). 
Although the current sample is too small to draw a conclusion, such a trend is to 
be expected if many of these MSPs are near accretion spin equilibrium, because 
pulsars with stronger magnetic fields should produce stronger oscillations and 
have lower equilibrium spin rates, other things being equal (see § 3). 

Classic kilohertz QPO pairs have been discovered in two accretion-powered 
MSPs, XTE J1807.4-294 (Markwardt et al. 2003a) and SAX J1808.4-3658 (Wi- 
jnands et al. 2003; see Fig. 3). The frequency separation of the QPO pair is con- 
sistent with the spin frequency of XTE J1807.4 but with half the spin frequency 
of SAX J 1808.4. The kilohertz QPO separation is consistent with the spin fre- 
quency or half of it in all stars in which burst oscillations have been detected (see 
Lamb & Miller 2004; Lamb 2005). These discoveries demonstrate conclusively 
that some kilohertz QPO sources have dynamically important magnetic fields 
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Figure 3. Power density spectrum of the X-ray brightness variations 
of the accretion-powered MSP SAX J1808.4-3658 on 18 October 2002, 
showing the 401 Hz periodic oscillations ("pulsations") at the star's 
spin frequency, the lower kilohertz QPO at 499ib4 Hz and the upper 
kilohertz QPO at 694±4 Hz (Wijnands et al. 2003). 



and that the spin of the star plays a central role in generating the QPO pair, 
confirming two important predictions of the sonic-point beat-frequency model 
(MLP98). But they also show that the original model is incorrect or at least 
incomplete, because it cannot explain a frequency separation equal to half the 
spin frequency (Lamb &: Miller 2004). 

A modified version of the sonic-point beat-frequency model (Lamb & Miller 
2004; Lamb 2005) attributes the lower kilohertz QPO to interaction of X-rays 
from the stellar surface with vertical motions of gas in the disk excited by the 
star's radiation or magnetic field at the radius where its spin frequency resonates 
with the vertical epicyclic frequency. The lowest-order linear resonance of this 
type would generate vertical oscillations with half the spin frequency and could 
produce QPO frequency separations equal to the star's spin frequency or half of 
it (Lamb & Miller 2004; Lamb 2005). 

The new results exclude models (see Lamb 2003) in which the frequencies 
of the kilohertz QPOs are various relativistic precession frequencies because 
spin frequencies several times higher than those observed would be required 
and because such models cannot explain the commensurability of the frequency 
separation of the QPOs with the stellar spin frequency (Lamb 2003; Lamb 2005). 



3. Evolution of Neutron Star Spins and Production of Millisecond 
Rotation-Powered Pulsars in LMXBs 

Current spins of neutron stars in LMXBs. — Accretion from a disk will spin up 
a slowly-rotating neutron star on the spin-relaxation timescale (Ghosh &; Lamb 
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1979, hereafter GL79; Ghosh & Lamb 1992, hereafter GL92; Lamb & Yu 2004) 

^ 2..,,.//[M(0Mr„)'/^l ~ m'yr (^) , (1) 

where fspinj M, and / arc the star's spin, mass, and moment of inertia, M is 
the accretion rate onto the star (not the mass transfer rate), is the angular 
momentum couphng radius, a is 0.23 if the inner disk is radiation-pressure- 
dominated (RPD) or 0.38 if it is gas-pressure-dominated (GPD), and in the last 
expression on the right the weak dependence of tspin 

on M, I, and the star's 

magnetic field has been neglected. 

The current spin rates of neutron stars in LMXBs reflect the average ac- 
cretion torque over a period ~ igpin 

or longer. Determining this average torque 

is complicated by the fact that the accretion rates and magnetic fields of these 
stars vary with time by large factors and that the torque can decrease as well 
as increase the spin rate. While a few neutron stars in LMXBs accrete steadily 
at rates ~ Me, most accrete at rates ~ 10~^-10~^M£; (Hasinger & van der 
Khs 1989; Lamb 1989; van den Heuvel 1992; MLP98) and many accrete only 
episodically (van den Heuvel 1992; Ritter Sz King 2001). The recently-discovered 
accretion-powered MSPs in LMXBs have outbursts only every few years, dur- 
ing which their accretion rates rise to ~ W~^Me for a few weeks before falling 
again to ^ 10"^Me (see Chakrabarty et al. 2003; Strohmayer et al. 2003). Also, 
there is strong evidence that the magnetic fields of neutron stars in LMXBs de- 
crease by factors ~ 10^-10^ during accretion, perhaps on timescales as short as 
hundreds of years (see Shibazaki et al. 1999; Bhattacharya &; Srinivasan 1995). 

If a star's magnetic field and accretion rate are constant, accretion will 
spin it up on a timescale ~ ispin 

to its equilibrium spin frequency u^^. This 
frequency depends on M, the strength and structure of the star's magnetic 
field, the thermal structure of the disk at r„j, and M (GL79; White & Stella 
1987; GL92). If a star's magnetic field and accretion rate change on timescales 
longer than ^spin, the spin frequency will approach i/gq and track it as it changes. 
If instead M varies on timescales shorter than tspm, the spin rate will fluctuate 
about the appropriate average value of Ucq (see Eisner, Ghosh, & Lamb 1980). 
Thus i^eq and its dependence on B and M provide a framework for analyzing 
the evolution of the spins of neutron stars in LMXBs. 

Figure 4 shows for five accretion rates and dipole magnetic fields = 
3.2 X 10^9 (PP) G ranging from 10^ G to 10^^ G. The lines are actually 
bandsj due to systematic uncertainties in the models. The lines for M = Me 
and M = O.lMg have jumps where the structure of the disk at the angular 
momentum coupling radius changes from RPD (lower left) to GPD (upper 
right); in reality the transition is smooth. For M ^ O.OIM^', the disk is GPD 
at rm even if the star's magnetic field is 3 x 10^ G. The effects of the stellar 
surface and the innermost stable circular orbit (Lamb &; Yu 2004) are not shown. 

The properties of the 16 known MSPs in LMXBs (Table 1) arc consistent 
with spin-up by accretion to spin equilibrium if their magnetic fields are between 
3 X 10'' G and 3 x 10^ G and their time-averaged accretion rates are between 3 x 
W^^Me and Me- The spin rates and visible pulsations of the accretion-powered 
MSPs are understandable if they have magnetic fields ~ 3 x 10^ G and have been 



8 



Lamb k, Yu 



T 



T 



10-10 - 



®Binary 



10-12 




> 



- 10-14 



> 




o 



10-18 - 



■rH 

CD 



10-20 



10-22 



0.001 0.010 0.100 1.000 10.000 
Period (s) 

Figure 4. Spin-evolution diagram. Lines sloping downward to the 
right show the P-P relation for magnetic dipole braking by a field 
with the strength indicated. Lines sloping upward to the right show 
the equilibrium spin period of a neutron star with the accretion rate 
indicated by the labels and a dipole field of the strength indicated by 
the downward-sloping lines. The dashed line sloping upward to the 
right shows where stars with a spin-down time equal to 15 Gy would 
lie. Data points are known rotation-powered pulsars; those of pulsars in 
binary systems are encircled. Data from Hobbs & Manchester (2004). 

spun up to accretion spin equilibrium by accretion at rates ~ \^~'^Me- The 
higher spin rates of the other MSPs are understandable if they have magnetic 
fields <^ 10^ G and time-averaged accretion rates ~ lO-^Mg. The absence of 
MSPs with spin frequencies > 750 Hz is consistent with spin-up to accretion 
spin equilibrium if these stars have magnetic fields 3 x lO'^ G and average 

accretion rates <^ 10"^ M^:. These fields and rates are consistent with the other 
observed properties of individual neutron stars in LMXBs (MLP98; Psaltis & 
Chakrabarty 1999; Chakrabarty et al. 2003). Alternatively, these stars may 
have stopped accreting before reaching spin equilibrium or been spun down as 
accretion ended. 

Based on the limited information then available, some authors (Bildsten 
1998; Ushomirsky, Cutler, & Bildsten 2000) speculated that neutron stars in 
LMXBs have negligible magnetic fields and spin frequencies in a narrow range, 
with many within 20% of 300 Hz. Such a distribution would be difficult to ex- 
plain by accretion torques and was taken as evidence that gravitational radiation 
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plays an important role. We now know (sec § 2) that most if not all neutron 
stars in LMXBs have dynamically important magnetic fields, that the observed 
spins of neutron stars in LMXBs are approximately uniformly distributed from 
< 200 Hz to > 600 Hz, and that production of gravitational radiation by uneven 
heating of the crust or excitation of r-waves is not as easy as was originally 
thought (Ushomirsky et al. 2000; Lindblom Sz Owen 2002). The spins of neu- 
tron stars in LMXBs may be affected by gravitational radiation, if their magnetic 
fields are weak enough, but there is at present no evidence for this. 

Production of millisecond rotation-powered pulsars. — The initial proposals 
that neutron stars are spun up to millisecond periods by accretion in LMXBs 
assumed that they accrete at rates ~ Me throughout their accretion phase (see 
Bhattacharya &i van den Heuvel 1991) and implicitly that accretion then ends 
suddenly; otherwise the stars would track to low spin rates as accretion ends. 
This simplified picture is sometimes still used (see, e.g., Arzoumanian, Cordes, 
&; Wasserman 1999), but — as noted above — most neutron stars in LMXBs ac- 
crete at rates <C Me and many accrete only episodically. The real situation is 
therefore more complex. 

The initial spins of rotation-powered MSPs recycled in LMXBs are the 
spins of their progenitors when they stopped accreting. These spins depend 
sensitively on the magnetic fields and the appropriately averaged accretion rates 
of the progenitors when accretion ends. Comparison of the equilibrium spin- 
period curves for a range of accretion rates with the P-P distribution of known 
rotation-powered MSPs (Fig. 4) suggests four important conclusions: 

(1) The overall P P distribution is consistent with the distribution expected 
for spin-up to accretion spin-equilibrium and generally supports the models of 
disk accretion by weak-field neutron stars used. In particular, the observed P- 
P distribution is consistent with the predicted distribution only if the accretion 
torque vanishes at a spin frequency close to the predicted fgq- 

(2) The accretion spin-equilibrium hypothesis predicts that MSPs should 
not be found above the spin-equilibrium line for M = Me and Uc = I, because 
this is a bounding case. The observed P-P distribution is consistent with the 
RPD model of the inner disk that was used for M ^ O.IMe, except for two 
pulsars discovered very recently in globular clusters (B1821— 24 [lower-left] and 
B1820-30A [upper-right]; Hobbs et al. 2004). Either the intrinsic P's of these 
pulsars are lower than shown or the RPD model of the inner disk does not 
accurately describe the accretion flow that spun up these stars. 

(3) The accretion spin-equilibrium hypothesis predicts that MSPs should be 
rare or absent below the spin-equilibrium line for M = 10~^ Me, because stars 
accreting at such low rates generally will not achieve spin equilibrium during 
their accretion phase. The observed P-P distribution is consistent with this 
prediction. 

(4) The MSPs near the 15 Gyr spin-down line were produced in situ by 
final accretion rates ^ 3 x 10~^ Me rather than by spin-up to shorter periods by 

accretion at rates ^ 3 x 10~^ Me followed by magnetic braking, because braking 
would take too long. This result accords with the expectation (sec above) that 
most neutron stars in LMXBs accrete at rates -C Me toward the end of their 
accretion phase. 
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4. Concluding Remarks 

The RXTE mission has discovered 5 accretion-powered and 13 nuclear-powered 
MSPs, with frequencies ranging from 185 Hz to 619 Hz. Nine of these MSPs 
produce kilohertz QPOs, and the kilohertz QPOs of another dozen neutron 
stars in LMXBs indicate that they also have 200-600 Hz spin frequencies. The 
RXTE observations indicate that the MSPs in LMXBs have magnetic fields 
ranging from 3 x 10^ G to 3 x 10^ G. If they still have magnetic fields this strong 
and 200-600 Hz spin rates when accretion ends, these neutron stars are likely to 
become radio-emitting MSPs. These discoveries strongly support the hypothesis 
that neutron stars in LMXBs are the progenitors of rotation-powered MSPs. 

The current spin rates of neutron stars in LMXBs reflect the average of the 
accretion torque over the current spin-relaxation timescale, which typically is 
^ 10^ yr. It is difficult to compute the expected average accretion torque over 
this period, because the torque-averaged accretion rate is uncertain. If accretion 
torque theory can be validated, measurements of current spin rates can provide 
estimates of the average accretion rate over this period to compare with binary 
evolution calculations. 

The properties of the 16 known accretion- and nuclear-powered MSPs in 
LMXBs are consistent with spin-up by accretion to spin equilibrium if their 
magnetic fields are between 3 x 10^ G and 3 x 10® G and their time-averaged 
accretion rates are between 3 x 10~^Me and Me- The spins of neutron stars in 
LMXBs may be affected by gravitational radiation if their magnetic fields are 
weak enough, but there is at present no evidence for this. 

The initial spins of recycled rotation-powered MSPs reflect their magnetic 
fields and accretion rates when accretion ended. The P P distribution of rotation- 
powered millisecond pulsars indicates that their initial spins are set by spin-up 
to accretion spin equilibrium. Many are well below the spin-equilibrium line 
for accretion rates Id^'^M^ but have very long magnetic-braking spin-down 
times, indicating that they had accretion rates ^ Me when accretion ended or 
never reached spin equilibrium. 

After more than two decades of effort, accreting MSPs in LMXBs have at 
last been found. The discovery of these stars promises important advances in 
our understanding of the evolution of neutron stars in LMXBs and the formation 
of radio-emitting, rotation-powered MSPs. 

We thank L. Bildsten, D. Chakrabarty, P. Kaaret, M. van der Klis, M. C. 
Miller, D. Psaltis, and J. Swank for helpful discussions. This research was sup- 
ported in part by NASA grants NAG 5-12030 and NAG 5-8740, NSF grant 
AST 0098399, and the funds of the Fortner Endowed Chair at the University of 
Illinois. 
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